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PREFACE 


This research was conducted according to the principles enunciated in the 
“Guide for Laboratory Animal Facilities and Care," prepared by the National 
Academy of Sciences, National Research Council. 
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SECTION 1 


INTRODUCTION 


Since the mid 60's, extensive studies have been carried out by the Armed 

Forces Radiobfological Research Institute (AFRRI) and other military organizations 

to examine the effects of intermediate to high doses of ionizing radiation on 

mammals. The intent of these studies has been to acquire information on lower 

order mammals, which might be extrapolated to man. Particularly extensive studies 

| have been carried out using the rhesus monkey (Macaca mulatta) and the AFRRI TRIGA 

reactor exposure facility. The purpose of the work described in this report is to 

duplicate these studies py analytical means. This is done in order to provide 

spectral detail not available from common integral measurement techniques and to 
obtain dosimetry data for locations not am2nable to measurements. 
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‘ Studies of effects of ionizing radiation on the rhesus monkey have involved 
post-exposure performance testing and assessments of mortality incidence. These 
effects have been correlated with dosimetric quantities which are thought to be 


: relevant to the effect in question if not to have a causai relationstiip. Dose to wenn 
the mid-head has generally been associated with degradation in performance or with ay 


incapacitation, while dose to a point at the mid-thorax has generally been 
associated with mortality incidence. 
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It ts generally recognized that the threshold dose which has a causal 
relationship to mammalian mortality is that which is deposited in the active 
i marrow. |!) It has also been theorized that incapacitation may be related to some 
is dose other than mid-head. Of course the reason why the mid-head and mid-thorax 
, doses have remained in vogue for correlation purposes is that they represent well 


defined and relatively unambiguous locations at which all types of measurements 
can be conveniently made. On the other hand it would require a distribution of 
detectors throughout an accurate phantom or cadaver to measure dose to marrow. NS 
Such experiments are difficult and inconvenient to perform and are restrictive as ~ 
to the types of detectors that can be used, due to size limitations. 
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Because of the practical limitations on dose measurements, analytical methods 
for obtaining the dosimetric quantities of interest have been examined with some 
success. In an earlier study'?) neutron and gamma spectra were obtained in simple 
geometries by experimental and analytical means and compared. In the study 
described herein this work has been extended to include analysis of radiation 
transmission within a detailed rhesus monkey phantom. 


The phantom has been developed based on anatomical measurements of a rhesus 
monkey cadaver. These measurements have been incorporated into a three-dimen- 
sional mathematical model through the use of combinatorial geometry. This is a 
system by which complex objects may be constructed analytically through the 
combination of simple shapes. Combinatorial geometry is interpreted for radiation 
transport analysis using the three-dimensional Monte Carlo code morse. (3) 


Calculations of the transport of neutron and gamma radiation within the 
monkey phantom due to such radiation externally incident on the phantom, 
particularly that of the AFRRI TRIGA reactor, have been performed and are 
described herein. The calculations include the mid-head, mid-thorax and active 
marrow aS analytical detector locations with exposure in two TRIGA exposure room 
configurations. A companion set of calculations for the mid-head and mid-thorax 
locations has been performed for a simple phantom. This has been done in order to 


allow validation and verification of the analytical method by means of comparison 
with experiment. 


The balance of this report describes the development of the realistic 
(complex) rhesus monkey phantom and its simple analogue. It also describes the 
calculations of radiation transport in the complex and simple phantom and presents 
the results in graphical and tabular form. Finally, the report describes the 


results of validation and verification studies of these calculated values as 
performed to date. 


SECTION 2 
RHESUS MONKEY PHANTOM 


This section of the report describes the development of specifications for a 
complex three-dimensional rhesus monkey phantom and its simple analogue. 


2.1 Complex Geometry Phantom 


Science Applications, Inc., has developed a complex computational model of an 


adult rhesus monkey in three dimensions, including detailed representations of 


lung, skeletal and soft tissue. The model has been produced using the 
combinatorial geometry module of the MORSE Monte Carlo radiation transport code 
and is represented pictorially herein using the GIFT code. (4) Combinatorial 
geometry is a system for representing complex shapes by means of combinations of 
simple shapes or "bodies" understood individually by the computer program. The 
bodies are assembled into regions, such as tne Jungs, skeleton and soft tissue of 
an anthropomorphic phantom. Certain limitations apply to this process which 
should be explained to the reader before describing the model and its basis. 
First, the process of computing radiation transport through the model media 
requires that the crossing of all boundaries between bodies be accounted for. 
This must be done in order for the code to know the location of the radiation 
particle and also to know the nature of the material through which it has passed. 
Such an accounting process is time consuming. Thus, in order for the transport 
program to operate efficiently the number of bodies which are used to represent an 
object should be kept to a minimum consistent with a reasonable level of physical 
accuracy. Second, the shapes of the bodies understood in combinatorial geometry 
have some limitations. The most important of these are 1) that an ellipsoid must 
be rotational about one axis and 2) that elliptical cylinders and truncated 
elliptical cones must have their upper and lower bases in proportion to one 
another, i.e., the ratio of the minor to major axis must be the same for the upper 
and lower bases. With these limitations in mind the monkey phantom model can now 
be described. 


ee Be eT Re i Ca al 
* 


a 
MO Ne De ae ee oe gat 


’ 
la 


ALL 


id 


2.1.1 Phantom Dimensions 


The monkey phantom has been developed based on anatomical measurements of a 
rhesus monkey cadaver. Such measurements were made from true scale photographs 
portraying the upper surfaces of lateral cross sections of the monkey, taken at 
approximately 3 cm intervals. Dimensions of the monkey, taken before and after 
sectioning, are given in Tables 1 and 2, along with equivalent values from the 
phantom model. 


The model was designed to correspond to the plan view of the photographed 
sections and the sum of their thicknesses. Consequently, given the effect of 
roundoff error and irregular sections, the model is two centimeters shorter, crown 
to tailbone, than the measurement of the intact subject. The lengths of the limbs 
were modeled to correspond to internal measurements, joint to joint, and are 
similar to but not exactly the same as the limb lengths measured externally. 
Hands and feet have been omitted from the model, since they represent little mass 
and contain little material of dosimetric interest. The partial body lengths are 
approximate, particularly since the model has no explicitly designated chin or 
mandible. Instead, the base of the chin in the model has been taken to be the top 
of the neck, and the crown-mandible length is taken to the maximum minor diameter 
of the elliptical head of the model. 


The model is slightly smaller than the prototype subject based on the 
external dimensions, including circumferences of the torso at two points. 
However, comparisons of the principal dimensions of the individual sections with 
those of the model indicate very close agreement between the two. The locations 
of these sections are depicted in Figures 1 (frontal view, Azimuth 0.0°, Elevation 
0.0°)) and 2 (side view, Azimuth 90°, Elevation 0.0°). They are described in 
terms of their height above the pad or the rump of the monkey, and are given as 
width (distance side to side) and depth (mid-sagital) dimensions normal to the 
length of the animal. 


Views of the cross sectional surfaces, looking downward from the level of the 
section (Azimuth 0.0°, Elevation 90.0°), are provided in Figures 3 through 19, 
which correspond to each of the sections listed in Table 2. It can be seen from 
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the Figures that the torso of the model is made up of concentric cylinders and 
cones, some elliptical. The shoulder region is constructed using additional 
elliptical cones which overlap those which make up the torso. 


Figures 20 through 23 show various views of the phantom. The phantom posture 
and orientation of the limbs have been designed to correspond to those of the 
monkey in a restraint chair. 


Figures 3 through 19 also show skeletal and lung detail corresponding to that 
found in the section photographs. The three dimensional representation of this 
detail is shown in Figures 24 through 27. 


The volumes of the lungs, skeletal and soft tissue regions are 124, 370, and 
4768 cubic centimeters, respectively, as determined by GIFT. Assuming that 
densities of equivalent human regions apply to the monkey (lung: 0.296, skeleton: 
1.40 and soft tissue: 1.04 g/cc'®)) the total mass of the model is 5.513 kg. This 
mass is consistent with rhesus monkeys of the general dimensions of the study 
subject as reported in published data. (6) According to these data, a monkey of 36 
to 42 months in age, having length (head/buttocks) of 47.4 cm, may be expected to 
have a mass of 5.2 kilograms. The subject monkey is between 51 and 53 cm in 
length, having other dimensions on a similar scale relative to that published and, 
thus, may be expected to have a slightly larger mass, as is indeed the case. 


Although the 5.5 kilogram mass can be justified as representative of the 
particular monkey being modeled based on available cadaver data, it is 
nevertheless at the high end for experimental subjects used at AFRRI as shown in 
Table 3. Therefore, SAI was directed by the sponsor to adjust the mass of the 
monkey to be approximately 3.4 kg,which is near the middle of the size range of 
subjects previously used at AFRRI and more of typical of modern subjects. After 
studying data on the rhesus monkey anatomy, as referred to previously, it was 
confirmed that the linear dimensions of the monkey were proportional to the cube 
root of the mass. The only complication in this was the fact that the 
proportionality for the head was found to differ from that of the remainder of the 
monkey in that head displayed less size variation with monkey mass. Constants of 
proportionality were determined for both monkey regions based on the two monkey 
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Table 1. Specifications of subject, as measured prior to 
sectioning, and model 


Measure (cm) 


2 Description Subject Model sgt 
i : : 
‘ Top of Crown - Tailbone 53 51 = eaey es 
Tip of Finger - Point of Elbow 23 14 (ulna) \ 
Elbow - Humeral/Scapula Joint 16 13 (humerus) victors... 

Point of Ankle - Knee 17 15 (tibia) ; e 
& Knee - Pelvis 13 14 (femur) ree : 
' Base of Chin - Mid Pelvis 33 32 ee 
K . a i™ 
# Crown - Lower Edge of Mandible 11 9 (head minor dia.) es 
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eae aah! 

Mass (Kg) 4.0 5.5 ic nae 

anes 

‘ ates 
iene 

nase 


sa 


~ 


Sgn ee 
2 Tie oe ce IT 


*relative to pad (seat) 


a eis 
5 ree 
‘ NC 
* wt: a ee *, 
¥ . say 1 
4 oe 
ti 
4 
4 
‘ 
: 12 
4 


[he eto eh Fo Sm Fe! . poe ane See ae Ee ee ie ae 
CART Rat bee 6 Pe Oho ooh ATO ORERES Cy Se 
OR NS ge ee Pe RECURRENCE R ORT ee ae 


Table 2. Dimensions of section upper surfaces 


Dimensions (cm) 


Subject Model 
Height above pad (cm) Width - Depth Width  - Depth 
3 9.1 5.3 91 6.8 
6 11.3 8.5 11.5 8.5 
9 12.0 10.0 11.5 9.5 
; 12 10.8 9.6 11.5 9.5 
15 10.5 10.6 10.3 10.3 
i 18 11.2 11.4 11.0 11.0 
F 21 11.9 11.8 11.8 11.8 
is 24 11.5 11.5 11.6 11.6 
27 11.1 11.5 11.4 11.4 
x 30 11.8 11.4 11.2 11.2 
: 33 @ mid depth 11.5 10.9 11.0 11.0 
i 33 max 15.2 13.4 
| 37 @ mid depth 12.1 B28) 98 8.1 
\ 37 max 13.2 12.8 
40 @ mid depth 8.0 5.5 8.4 6. 
5 40 max 9.3 9.0 
42 no data available 4.5 4.5 
P 45 8.0 11.5 7.9 9.8 
‘i 48 8.5 10.0 7.9 9.8 
i 50 5.0 6.5 5.4 5.6 
: Total Section Thickness 51 51 
y, 
. 
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Table 3. 


Specifications of age and mass for rhesus monkeys 


used in AFRRI experiments, 1966 to 1973 


' Report No. No. of Subjects 
SR66-2 131 
SP66-23 190 
SR68-16 6 
SR69-1 6 
SR69-2 6 
SR69-8 6 
; SR69-9 6 
SR69-14 6 
SR69-18 30 
SR69-21 14 
’ SR70-7 7 


SR70-11 
SR71-3 
SR71-6 


} $R71-10 6 
i SR72-14 8 
SR72-19 85 
SR73-1 88 
4 Subject-weighted average mass range: 
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Figure 1. Locations of cross sections, frontal view 
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Figure 3. Monkey phantom cross section, Z = 3.0 cm 
azimuth 020° elevation 90,0° 
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Figure 6. Monkey phantom cross section, Z = 12.0 cm 
azimuth 0.0° elevation 90.0° 
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Figure 7. Monkey phantom cross section, Z = 15.0 cm 
azimuth 0.0° elevation 90.0° 
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Figure 8. Monkey phantom cross section, Z = 18.0 cm 
azimuth 0.0° elevation 90.0° 


22 


eee ef 


“ 


= 4 % 
Oey 


. 


. 


ore 


Se s,s ee ee Te OE Se Te 2 os ~» dg he 
Waithesiy, 
aks, 


Baits, " 
Gn) At 


2 
tt 
a 
* 
G4 


5 
aa dy 


ht 
x 


Pale 


tate 


‘ Ayr 


| Soft Tissue | ‘ 
; his eee 


Skeletal Tissue 


Lung Tissue 


2 Figure 9. Monkey phantom cross section, Z = 21.0 cm 
8 azimuth 0.0° elevation 90.0° 
4 


Pe 23 


Pid td Ae 


. ‘ 
"v « N, <> 


te Se Om Me Mw 


24.0 cm 


> 
ize 


a 
cc 
fo) 
.- 
Pe) 
U 
3) 
n 
"” 
no 
id (=) 
wo ° 
3 vo Oo 
7) “ ov n 
3 7) 3 S = 
2) bane” " 
“ -— 7.) - oO Tt 
_ Yy - cS we N 
-— - | ol cc +» 
wo Yj = ie] 
» ~ 7. f=2) Qa > 
Ge w c wv 
Oo -_ 3 > — 
Wn wy a yu wo 
x x 
(Fa) Cc 
[9] 
= oO 
[o) 
° 
- © 
= 
= <= 
wo 
=. 3s 
eS 
Om oe 
- N 
uw 6 


Soft Tissue 


Skeletal Tissue 


NN 


Lung Tissue 


Figure 11. Monkey phantom cross section, Z = 27.0 cm 
azimuth 0.0° elevation 90.0° 
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Figure 12. Monkey phantom cross section, Z = 30.0 cm 
azimuth 0.0° elevation 90.0° 
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Figure 13. Monkey phantom cross section, Z = 33.0 cm 
azimuth 0.0° elevation 90.0° 
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Figure 16. Monkey phantom cross section, Z = 42.0 cm ees 
azimuth 0.0° elevation 90.0° s 
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Figure 17. Monkey phantom cross section, Z = 45.0 cm 
azimuth 0.0° elevation 90.0° 
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Figure 22. Monkey phantom, side view 
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sizes discussed in the reference, 3.5 and 5.2 kg, (6) and used to scale the 5.5 kg wpe 
phantom to the smaller size. The ratios of new to old phantom masses were 0.76 SENS 
and 0.61 for the head and for the remainder of the phantom, respectively. Tee 
Likewise, the location of the head midpoint relative to the base of the trunk was | 
adjusted from 46.5 cm to 39.7 cm to remove any gap resulting from the scaling 
procedure. In its final adjusted form the monkey phantom has a mass of 3.43 kg 
and has a head which is slightly larger, relative to the trunk or torso, than that 
of the original phantom. 
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RS 2.1.2 Elemental Composition 
¥ 
; A specification of elemental constituents of monkey tissues is not available. 
ey Therefore, the latest constituents determined by kerr?) to apply to human lung, 
3 skeletal and soft tissues have been used to represent the equivalent monkey 
a tissues until such time as alternative values become available. These 
constituents are listed as weight percents in Table 4. 

Bat 
: 2.1.3 Marrow Distribution 
m 

The active marrow distribution in the rhesus monkey was measured by S.T. 
Be Taketa, et al., and published in 1970. (8) The results of these measurements along 
- with those for other species are given in Table 5, which is reproduced from that 
RG report. Figure 28 shows this distribution pictorially, based on scintillation ee, 
- scanning of 99m. in an eviscerated monkey. This marrow distribution has been » t 
a used in the dosimetry calculations as shown in Figure 29, simplified to correspond ae 
a to the seven gross marrow regions given in Table 5. Some additional detail has see : 
a) been provided in the clavical-scapular region and in the legs increasing the total a 
. individual marrow regions to nine. The division of marrow between the clavical et 
es and scapulae are as shown in Table 5. Division of marrow between the upper femur tee 
ra and the remainder of the leg was accomplished by placing approximately half the Ne: 
7 femur marrow (6.7%) in the region near the pelvis and the remainder of the lower AS 
a limb marrow (13.3%) in the area of the knee joint as supported by the Taketa | 
“ reference. eae 
oe oe 
42 y 4 
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Table 4. Elemental composition for various components of reference man 
to be used as Rhesus monkey analogues 


5 Percent by Weight 
Element _tung_ Skeleton Soft Tissue 
H 10.21 7.28 10.51 
c 10.24 24.64 22.63 
q N 2.91 3.06 2.34 
‘ 0 75.63 46.88 63.69 
d Na 0.19 0.32 0.11 
Mg 0.007 0.11 0.013 
* p 0.080 5.03 0.13 
s 0.23 0.31 0.20 
i cl 0.27 0.14 0.14 
; K 0.20 0.15 0.20 
s Ca 0.009 12.07 0.024 
a Fe 0.037 0.008 0.006 
; 
} 
: 
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Figure 28. Scintillation scanning of dal in an 
eviscerated monkey 


Scintillation scanning of Pate in an "“eviscerated" monkey (A) and 


of disarticulated bones (B) from the same animal. Radioactivity was 
of sufficient intensity to darken film in 1, skull and maximilla 
(positioned on left side with the anterior portion facing bottom of 
page); 2, humeri; 6, femora; 7, tibiae; 10, scapule; 11, lumbar anc 
12, thoracic vertebrae; and 14, hip bones; but not in 3, radii; 4, 
ulnae; 5, mandible; 8 fibulae; 9, clavicles; 13, cervical vertebrae; 
15, feet and ankles; 16, hands and wrists; and 17, gacrococcygeal 
vertebrae. (Drawings were reconstructed from the 97™TC scans and 
radiographs of the skeleton.) 
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2.1.4 Monkey in Experimental Configuration 


A primary mode of exposure for the monkey is more or less centrally located 
within a plywood box, restrained in a chair made predominateiy of lucite plastic, 
with its back to the reactor. This disposition of the monkey is shown in Figure 
30. The general dimensions of the chair and box are given Figures 31 and 32, 
respectively. The elemental constituents and densities of these materials are 
given in Table 6, (2 The density of the box was determined by dividing the 
measured mass of the box (44.7 kg) (9) by the volume of its sides, as determined by 
external measurements. 


The box and chair sit on a platform (not modeled) on the floor of the 
exposure room. The monkey position in the chair is adjusted to place the 
mid-thorax location at the same height as the reactor core center. The mid-thorax 
location is generally taken to be 3 cm below the end of the sternum. However, in 
the case of the phantom the sternum is not modeled separately. Therefore, the 
location chosen was alternatively defined as being below both lungs and ribs and 
between 19 and 20 cm below the mid-head in a 3.4 kg monkey. In the scaled 3.4 kg 
phantom the location of the mid-thorax was taken to be 20 cm above the base, which 
is 19.66 cm below the mid-head point. This location satisfies the conditions of 
being below both lungs and ribs. 


For irradiations in exposure room one, the box and chair are positioned atop 
a wood stand to bring the mid-thorax to 120 cm above the floor. Thus, in both 
exposure rooms the mid-thorax is aligned with the height of the center of the 
reactor core. Depending on the exposure room in which the experiment is to be 
performed, the mid-thorax may be either 86 or 92 cm above the platform. The exact 
location chosen affects the relative location of the monkey within the box. For 
the purposes of this study the mid-thorax was taken to be located at 86 cm above 
the bottom of the box, as shown in Figure 33. 


2.2 Simple Geometry Phantom 


An exact physical mockup of the rhesus monkey phantom would be difficult if 
not impossible to construct for the purpose of experimentally validating the 
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Table 6. Elemental constituents of the 
lucite chair and plywood enclosure 


Element 


Lucite Density: 


Plywood Density: 


Constituent wt % 


Lucite P} ywood 
8.1 6.2 
60.0 44.5 
31.9 49.3 
1.2 g/cm? 
0.51 g/cm? 
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Figure 30. Chair and box with complex phantom 
) 
azimuth 175.0 elevation 5.0° 
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Orientation of complex phantom and chain within 


the plywood enclosure 


Figure 33, 


radiation transport calculations described in this report. Therefore, a simple 
geometry phantom was developed, which could easily be replicated with a physical 
model for use in experimental studies. The dimensions of the simple phantom were 
chosen to generally correspond to those of the unscaled phantom head, neck and 
trunk, minus the shoulders. The specifications arrived at in cooperation with the 
sponsor (who will actually build the phantom) are as follows: 


Trunk: cylinder 

height: 33.0 cm 

diameter: 10.2 cm 

composition: TE liquid in a lucite cylinder with 0.125" thick walls 
Neck: cylinder 

height: 4.5 cm 

diameter: 4.5 cm 

composition: A-150 plastic 
Head: sphere 

diameter: 9.5 cm, truncated to a height of 8.0 cm, with the flat 

side sitting on neck 
composition: A-150 plastic 


The constituents of lucite were given previously in Table 6. The constituents of 
TE (Tissue Equivalent) liquid and A-150 plastic are given in Table 7, (10) 


The simplified phantom was located within the restraining chair and exposure 
box as shown in Figure 34. The mid-head location of the simple phantom within the 
box was taken to be the same as that of the complex, with the mid-thorax specified 
to be 22.5 cm below. The exposure configuration thus obtained is shown in Figure 
35. 
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Table 7. Elemental constituents of the simplified phantom 
Constituent Wt% 
Element Tis. Equiv. Liquid A-150 Plastic 
H 10.2 10.1 
C 12.0 77.6 
N 3.6 3.5 
0 74.2 5.2 
Ca 1.8 
F 1.7 
Tissue Equivalent Liquid Density: 1.08 g/cm® 
A-150 Plastic Density: 1.11 g/cm? 
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Figure 35. 
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SECTION 3 
FLUENCE AND KERMA CALCULATIONS 


Calculations of neutron and gamma_ radiation transport within the 
phantom-chair-box system were performed using the MORSE Monte Carlo code in the 
adjoint mode. In that mode particles are started at the detector position within 
the phantom and are followed through an inverse of the applicable physical 
processes until they leak from the system and are tallied by energy and direction. 
This process is achieved by inverting the transport cross section matrix and then 
performing the transport calculation much as if it were in the normal forward 


mode. 


The value of the adjoint approach is that the results of the calculation are 
the equivalent of the detector response as perturbed in energy and angle by the 
surrounding media. Thus, if one can make the assumption that those media do not 
perturb a radiation field beyond their own confines, the results of the 
calculation may be used to determine fluence values (E’,9') internal to the 
system placed in any arbitrary location or orientation within the radiation field, 
according to the expression: 


OE", 2") =f ff 2 (S.E,0) *(5,E8", 2+ 9) Ne Md Qdeds 
S E Q 


where (S,E,2) is the free-field radiation fluence impinging on the surface at 
S, differential in energy (E) and angle (2), * is the adjoint fluence, born at 
the detector with energy E' and angle &', and leaking from the surface at S with 
energy E and angle 2. Finally, N is a unit vector normal to the surface at S. 


The above integral process is often referred to as forward-adjoint coupling 
and is accomplished using the VCS code system, (11) This code system is capable of 
associating each individual Monte Carlo adjoint leakage history with the 
appropriate free-field angular fluence obtained by one-, two or three-dimensional 
calculation methods. Because these leakage histories are fully correlated with 


% 


the description of each individual particle started at the detector site, the code 
system is capable of calculating either the energy-differential fluence at the 
detector, which results from the incident fluence, or the associated KERMA or any 
other appropriate quantity. 


Calculations of neutron and gamma ray transport were performed for detector 
locations at mid-head and mid-thorax for both the complex monkey phantom and its 
simplified analogue. Calculations were also performed for detectors corresponding 
to the marrow distribution in the complex phantom, as previously described. These 
calculations were performed in two parts, the first starting approximately 100,000 
neutrons and the second starting 100,000 gamma rays, each, for the mid-head and 
mid-thorax locations. This number of starting particles has been found to 
produce statistically reliable results (fractional standard deviations of 0.10 or 
lower in each of the three fluence components: neutrons, neutron-induced gamma 
rays and gamma rays, when each is KERMA-weighted and summed over energy and 
angle). For each of nine marrow regions 20,000 of each radiation type were 
started, giving a total of 180,000 starting particles of each type for the marrow 
as a whole. 


Cross sections were obtained from the DNA Few Group Cross Section Library 
(pLc-31). 12) This library contains coupled sets of multigroup transfer cross 
sections in 37 neutron and 21 gamma energy intervals as shown in Table 8. The 
neutron cross sections were collapsed into the few group format using an 
assumption of in-group fluence proportional to the inverse of the energy. The 
thermal group was Maxwellian (300° K) weighted, The gamma ray cross sections were 
collapsed using uniform in-group weighting. The angle of scattering was treated 
using a Ps Legendre polynomial expansion. 


The results of the calculations have been retained in their fully detailed 
form as previously described and are available for reuse at anytime. They are 
depicted here as adjoint leakage scalar fluence from a soft tissue KERMA detector. 
The KERMA values are provided in Table 9 and correspond to the components of ICRU 
muscle’ 10) , Neutron KERMA were calculated using the MACK code! 13) | while those 


8 GEES 14 +) Fe! 


for gamma rays were obtained from DLC-31. The adjoint KERMA-weighted leakage 
scalar fluence K*(E) is actually nothing more than the detector KERMA values 
perturbed by transport through the system and may be defined as: 


k*(E) SffSf K(E") @*(S,E + E', Q> Q')N-QdQd Q'dE'dS 
S E'Q2'2 

where K(E') is the KERMA value at detector energy E'. ®* represents the adjoint 

leakage fluence from an unweighted or uniform detector response used as the 

adjoint source, which is in turn weighted by K(E'). 


The six figures beginning with Figure 36 depict the adjoint neutron scalar 
fluence for neutron and neutron-induced gamma ray KERMA and the adjoint gamma ray 
scalar fluence for gamma ray KERMA, applicable to each phantom type and dosimeter 
location. These data are presented as group fluences and are similarly tabulated 
in the Appendix. This format is used to present the adjoint scalar fluence in a 
form equivalent to the KERMA, which are used as the adjoint source. This format 
is also used later in the report to present the transmitted fluence spectra, 
because it allows for simple convolution with group KERMA to obtain total KERMA 
values. The horizontal scale is linear above 1 MeV to provide better resolution 
for the reader in this dosimetrically important range. For the complex phantom 
the results for mid-head, mid-thorax and marrow locations are plotted together for 
each radiation component. The KERMA values are coplotted for comparison purposes, 
since, as mentioned previously, the adjoint fluences may be regarded as 
transport-weighted KERMA. For the neutron and gamma components the adjoint 
leakage fluences are larger for the mid~head than for the mid-thorax location, 
reflecting the relative masses of material surrounding the two sites. Because the 
marrow is distributed among shallow and deep sites the aggregate marrow response 
shows an irregular relationship with but is generally between those of the other 
two locations. The behavior of gamma ray production by neutrons as a function of 
neutron energy in the complex phantom is dominated by a minimum caused by a 
decrease in the likelihood that neutrons of increasing energy, but still below the 
inelastic threshold, will undergo the series of scattering events required for 
thermalization and capture without leaking from the system. Likewise, geometry 
plays a part in the production of gamma ray KERMA by incident neutrons. The 
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mid-thorax presents the most favorable geometry by surrounding the dosimetry point 
with production media in a nearly uniform manner. The mid-head location is a bit 
more toward an extremity and so receives fewer of the volume-generated gamma rays. 
That marrow component which is located at shallow locations receives even less, 
since the gamma rays can reach such marrow only from half-space. Hence, the 
aggregate marrow neutron-gamma response resembles that of the mid-head rather than 
the mid-thorax. 


The adjoint fluences from the simple phantom show characteristics similar to 
those of the complex phantom, although the neutron adjoint fluences are slightly 
lower for the simple phantom. The neutron-induced gamma ray responses are also 
similar. However, the minimum, which was so evident in the complex phantom 
response, is considerably less pronounced in the case of the simple phantom. The 
simple phantom is not more massive overall than the complex phantom. However, it 
is more compact, particularly in the trunk, i.e., it is larger in diameter and 
shorter, which results in less leakage of neutrons subsequent to their penetration 
of the system. 


It is recognized that, while of some academic interest, the adjoint leakage 
scalar fluence may in fact be of little relevance to the understanding of the 
dosimetry characteristics of a complex system, such as the monkey-chair-box. This 
is particularly true when the dosimeter location within such a system favors a 
particular exposure direction and the experimental radiation field is itself 
highly directional. As a first order attempt to understand the non-isotropi:- 
nature of the system, it is useful to examine the adjoint fluence integrated over 
Yess than 4n steradians. To this end, the three figures beginning with Figure 42 
depict the adjoint leakage fluence from the mid-thorax location in the realistic 
(complex) phantom, integrated separately over 47 steradians and over the 
hemisphere facing the reactor (front 21) and that facing away from the reactor 
(back 27). This case is considered to be qualitatively representative of all 
cases studied here. However, these data are tabulated for all calculated 
configurations in the Appendix. 
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Table 8. Neutron and gamma-ray energy boundaries for the 
37-21 coupled neutron-ganma library 


Neutron Group (eV) Gamma Group (eV) 
Group # Energy Lethargy Energy Eff. Avg. 
1 1.96+7 ~0.675 1.40+7 1.20+7 
2 1.69+7 ~0.525 1.00+7 9.00+6 
3 1.49+7 -0.400 8.00+6 7.50+6 
4 1.4247 -0.350 7.00+6 6.50+6 
! 5 1,38+7 ~0.325 6.00+6 5. 50+6 
6 1.28+7 -0.250 5.00+6 4.50+6 
| 7 1.22+7 ~0.200 4.00+6 3.50+6 
8 1.11+7 -0.100 3.00+6 2.75+6 
| 9 1.00+7 0.000 2.50+6 2.25+6 
| 10 9.0546 0.100 2.0046 1.7546 
| ll 8.1946 0.200 1.50+6 1.25+6 
: 12 7.41+6 0.300 1.00+6 8.5045 
| 13 6.38+6 0.450 7.00+5 5.7545 
| 14 4.97+6 0.700 4.5045 3.7545 
| 15 4.7246 0.750 3.00+5 2.2545 
| 16 4.0746 0.900 1.5045 1.2545 
17 3.01+6 1.200 1.0045 8.50+4 
18 2.39+6 1.433 7.00+4 5.57+4 
: 19 2.3146 1.467 4.50+4 3.7544 
| 20 1.83+6 1.700 3.00+4 2.50+4 
21 1.1146 2.200 2.00+4 1.50+4 
| 22 5.5045 2.900 1.00+4 
| 23 1.58+5 4.150 
24 1.1145 4.500 
25 5.25+4 5.250 
26 2.48+4 6.00 
27 2.19+4 6.125 re! 
28 1.03+4 6.875 o. 
29 3.3543 8.000 oe 
30 1,.23+3 9.000 
31 5.83+2 9.750 
32 1.01+2 11.500 
33 2.90+1 12.750 
34 1.0741 13.750 
35 3.06+0 15.000 
36 1.13+0 16.000 
37 E1421 17.000 
Lower Bound 1.00-5 27.631 
*Read as 1.96x107. 
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Figure 36. KERMA-weighted adjoint fluence complex phantom (n) 
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Figure 37, KERMA-weighted adjoint fluence complex phantom (n-g) 
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Figure 38. 
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Figure 39, KERMA-weighted adjoint fluence simple phantom (n) 
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Figure 40. KERMA-weighted adjoint fluence simple phantom (n-g) 
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The neutron and gamma ray adjoint fluences of the system, which includes the 
complex phantom, chair and box, have slightly larger values in the direction 
facing the reactor, particularly at lower incident energies. This corresponds to 
exposure of the back of the monkey phantom itself. The bias in favor of this 
direction is not large and, since the mid-thorax is centrally located front to 
back, is probably due to the shielding contribution of the box. The effect is 
likely to be enhanced somewhat for the marrow, which is somewhat biased toward the 
posterior of the phantom. 


The neutron-induced gamma ray adjoint fluence, shown in Figure 43, shows less 
preference for the frontal exposure except at very low neutron energies. It is 
expected that these are the most affected by the box surrounding the monkey. 


It should be noted that the adjoint fluence integrated over 4n or even 2m 
does not reflect any variation with polar angle, which is significant, as will be 
described later in this report. 


3.1 Free-Field Coupling 


It is the ultimate goal of these calculations to be able to predict fluence 
and KERMA values in instances which are not practical for experimental dosimetry. 
To do this using the forward-adjoint coupling technique, described previously, 
requires that an applicable free-field fluence be defined. This has been done in 
a previous program using one and three-dimensional calculational techniques. (2) 
The one-dimension results were found to give good representations of the 
free-field scalar fluence spectra, although the absolute normalization and, 
perhaps, the angular distribution were not correct due to simplifications in 
source geometry. Because monitor data are available for these experimental 
exposures, the absolute normalization is not a critical issue. On the other hand 
no confirmation of radiation fluence angular distribution calculated by one 
dimensional techniques has been accomplished using experimental or 
multi-dimensional calculations. However, such multi-dimensional calculations are 
beyond the scope of this study. Thus, one dimensional energy and 
angle-differential fields have been chosen to represent the AFRRI exposure room 
radiation environments. 
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AFRRI has two exposure rooms, characterized as ER1 and ER2 in the previous 
study. Forward-adjoint coupling has been performed for each of these rooms in 
separate configurations. The ER2 configuration is that of the bare room, with the 
dosimetry point 100 centimeters from the reactor core centerline, 92 centimeters 
from the floor. This corresponds to case 5 in the previous study. The ER1 
configuration includes a two inch thick lead slab between the reactor and the 
dosimetry point, which is 100 centimeters from the core centerline and 120 
centimeters above the floor. This corresponds to case 7 in the previous study. 
The neutron and gamma spectra ray spectra calculated for these two locations are 
depicted in Figures 45 and 46, respectively. These have been normalized to 
produce the following measured soft tissue free-in-air (FIA) KERMA rates (rad-FIA 


per Kw-min) :(9) 


ER1 with 2" Pb @ 100 cm 
neutron 7.80 
gamma ray 2.22 

ER2 bare @ 100 cm 
neutron 10.8 
gamma ray 19.8 


Beside the obvious reduction in gamma rays caused by the lead it should be noted 
that the lead also causes a substantial shift in the neutron spectrum toward lower 
energies in the region of 1 MeV and above. It also results in substantial 
depletion of thermal and epithermal neutrons. It should be noted that 
measurements made at other distances from the reactor! 14) show a reduction of 
free-field intensity of ten percent or more at locations offset from the reactor 
centerline, corresponding to the mid-head location. However, no free-field values 
for the mid-head location are available at the 100 cm distance. Therefore, for 
the purposes of this report al} data are normalized to the intensity measured 100 
cm from the reactor midline. 


The VCS system calculates the transmitted fluence at the detector from an 
externally incident fluence, taking into account its full energy-angle detail. 
However, before examining the calculated transmitted fluences it is useful to 
examine the character of such transmission on the basis of transmission factors 
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calculated by dividing the KERMA-weighted transmitted fluence by that of the 
free-field. These transmission factors are presented in Table 10 for the complex 
and simple phantoms in each exposure room, along with their corresponding 
fractional standard deviation (FSD). 


Looking first at the simple phantom, there appears to be little difference in 
a either room between transmission to the mid-head and to the mid-thorax, at least 
: for neutrons and gamma rays and within the statistical uncertainties in the 
calculations. In the case of the neutron-induced gamma ray transmission, obtained 


“ by dividing the gamma ray KERMA produced by the neutron interactions in the 
‘ phantom by the neutron free-field KERMA, the transmission to the mid-thorax is 
in considerably greater than that to the mid-head. This is due to geometry 


considerations, as previously discussed. 


Looking further at the differences in transmission between the two rooms it 
can be seen that the incident gamma ray spectra do not differ sufficiently to 
produce significantly different transmission factors in the two rooms. On the 
other hand, the neutron spectra are sufficiently different to produce something of 
a paradox between the neutron and neutron-gamma ray transmission factors for the 
f two rooms. The absence of the lead shield in room 2 allows a large number of 
epithermal and thermal neutrons to be incident on the phantom in that room. These 
count for very little in so far as neutron KERMA is concerned but are very 
efficient in producing gamma rays in the phantom. Thus, the neutron-gamma KERMA 


: per unit free-field neutron KERMA are largest in room 2. However, even through 
i this is evidence of a softer neutron spectrum in room 2, the neutron transmission 
7 in that room is higher than in room 1. This seeming paradox is explained by the 


fact that, in the energy region most heavily weighted by the KERMA, it is the 
neutron spectrum in room 2 which is the hardest or most energetic and, hence, the 
most penetrating. The lead in room 1 softens the spectrum in the KERMA-weighted 
' region by repeated scatterings, while absorbing the low energy neutrons which have 
| little KERMA value. 


There is little qualitative difference between the behavior of transmission 
for the simple and complex phantoms, though there seems to be some evidence for 
: consistently less gamma ray transmission to the complex mid-head. This may be due 
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Table 10. Calculated transmission factors for the 
complex and simple monkey phantoms 


Exposure Room 1, @ 100 cm from core centerline, with 2 inches of lead 
intervening. 


Transmission Factor (FSD) 


Location Component Complex Simple 
Mid Head n 0.635 (.03) 0.572 (.03) 
n-¥ 0.040 (.07) 0.058 (.06) 
0.909 (.03) 0.971 (.03) 
Mid Thorax n 0.593 (.03) 0.572 (.03) 
0.952 (.03) 0.925 (.03) 
Active Marrow n 0.685 (.03) 
n-¥ 0.048 (.05) 
7 0.919 (.03) 


Complex Phantom 


Marrow Region/Component n n-¥ Y 
Skul] 0.493 (.07) 0.025 (.18) 0.756 (.06) 
Clavicles 0.705 (.07) 0.036 (.13) 0.851 (.06) 
Scapulae 0.635 (.07) 0.042 (.13) 0.967 (.06) 
Upper Limbs 0.620 (.07) 0.038 (.14) 0.916 (.06) 
Ribs 0.701 (.07) 0.049 (.12) 0.903 (.06) 
Vertebrae 0.749 (.06) 0.063 (.11) 1.036 (.06) 
Pelvis 1.037 (.06) 0.057 (.14) 0.953 (.05) 
Upper Femur 0.987 (.06) 0.059 (.15) 0.961 (.06) 
Lower Limbs 0.299 (.09) 0.027 (.15) 0.711 (.06) 
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Table 10. Calculated transmission factors for the *e 
complex and simple monkey phantoms (continued) * 


ees 


Exposure Room 2, @ 100 cm from core centerline, bare room. 


Transmission Factor (FSD) 


# Location Component Complex Simple 
i Mid Head n 0.738 (.04) 0.661 (.04) 
Be n-7 0.069 (.06) 0.090 (.05) 
0.920 (.03) 0.968 (.03) 
Mid Thorax n 0.718 (.04) 0.695 (.04) 
n-¥ 0.136 (.05) 0.131 (.05) 
0.972 (.03) 0.957 (.03) aE 
Active Marrow n 0.722 (.03 = 2 
n- 0.087 (.05) oy ati 
Y 0.935 Ba 
we * 
SN 
MN es 
Soe 
Complex Phantom i ee 
Marrow Region/Component n n-¥ Y Stee 
» Skul] 0.541 (.09) 0.047 (.14) 0.728 (.06) Mt 
; Clavicles 0.859 (.08) 0.074 (.10) 0.887 (.06) Beat 
Scapulae 0.749 (.08) 0.093 (.12) 0.971 (.06) ae 
Upper Limbs 0.735 (.08) 0.074 (.12) 0.925 (.06) rae 
Ribs 0.793 (.08) 0.093 (.10 0.916 (.05) 
" Vertebrae 0.836 (.07) 0.107 (.10) 1.058 (.05) 
Pelvis 1.164 (.07) 0.106 (.11) 1.006 (.05) 


Upper Femur 
Ps Lower Limbs 


0.994 (.07) 
0.337 (.11) 


0.104 (.13) 
0.050 (.14) 


1.020 (.05) 
0.696 (.06) 


to the skeletal material of the skull, which has a higher density than soft 
tissue, thus affecting gamma ray transport, but little difference with soft tissue 
in hydrogen density, thus leaving neutron transport relatively unaffected. In 
fact the elliptical geometry of the complex head appears to provide for increased 
transmission of neutron KERMA with an attendant decrease in neutron-induced gamma 
ray production, compared to the simple phantom. 


The mean value of the transmission of KERMA to marrow locations in the 
complex phantom is similar to that for the mid-head, though the neutron 
transmission is consistently larger. However, transmission values for the various 
marrow regions show that the mean value is made up of widely differing components. 
Those range from the low transmission factors for the best shielded components, 
such as the lower limbs and the skull, to transmission factors of the order of 
unity for the least well shielded locations, i.e., the vertebrae, pelvis and upper 
femur. Buildup in the phantom and scatter from the surrounding chair and box play 
a role in producing such large transmission factors for marrow locations closest 
to the reactor. The low transmission to the lower limbs is also understandable. 
However, the low transmission to the skull, particularly as compared to another 
location similar in depth and distribution, the ribs, is subject to question. it 
is possible that the box is responsible. However, if this is so then a premium is 
placed on the correct angular distribution incident on this location. As noted 
previously, the free-field angular distribution used in these calculations has not 
been verified. 


Figures 47 through 62 depict the transmitted neutron and gamma ray spectra at 
the various dosimetry points in the complex and simple phantoms in each exposure 
room. Data appearing in these figures are tabulated in the Appendix along with 
that for all detailed marrow locations. In each figure the free-field (FF) scalar 
fluence is depicted as a thick line. The mid-head fluence is represented by a 
thin solid line, the mid-thorax fluence by a thin dotted line and, in the case of 
the complex phantom, the marrow fluence is represented by a thick dotted line. 


Figures 47 through 50 depict the fluences in the complex phantom in exposure 
room one. The transmitted neutron fluences in Figure 47 show substantial 
reduction below the free-field above one MeV. Below that energy the phantom tends 
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to increase the transmitted fluence relative to the free-field, particularly in 
the thermal region. The neutron-induced gamma fluence shown in Figure 48 is 
noteworthy mainly for the large peak between 2 and 2.5 MeV, which is the energy 
range of the hydrogen capture gamma ray. In this one group the fluence from the 
neutron-induced gamma rays approches and even exceeds that of the room free-field 
gamma rays as shown by the thick solid line. The gamma ray transmission shown in 
Figure 49 is noteworthy for the large buildup of gamma ray fluence below one MeV. 
This buildup is particularly noticeable for the mid-thorax location, which is 
probably due to the fact that, while the torso and head of the complex phantom are 
of similar diameter, the skull material tends to screen out low energy gamma rays, 
as does the nearby box. 


SVAN MECC MS 


Figure 50 shows the total gamma ray fluence in the complex phantom co-plotted 
with the free-field. The large enhancement of the gamma ray fluence in the range 
about 100 keV and 1 MeV are very noticeable. At higher energies there is some 
enhancement due to the lack of attenuation by the phantom for incident gamma rays 
and some production of high energy gamma rays by incident neutrons, particularly 
in hydrogen and nitrogen. The marrow fluence appears to be more in line with that 
of the mid-head than the mid-thorax, tirough there is considerable statistical 
uncertainty in the fluence in individual energy groups. 


Figures 51 through 54 provide information for the complex phantom in exposure 
room two similar to that described above for room one. The general 
characteristics of the transmitted fluences in the two rooms are similar even 
though differences exist in the free-field spectra. Thus, the build up of the low 
energy neutrons and gamma rays in the phantom is less noticeable, since the 


PMS | = =A 


' 
aL Ae 


r 


om 
Ae 


: free-field contains more to begin with. Further, the effect of the 
- neutron-induced gamma rays in the phantom is less pronounced because of the low 
i ftree-field neutron-gamma ray ratio. 


Figures 55 through 58 depict the transmitted fluence calculated via VCS for 
the simple phantom in exposure room one. In Figures 55 and 57 it can be seen that 
for the high energy (E > 1 MeV) portions of the transmitted neutron and gamma ray 
spectra there is no consistent difference between those for the mid-head and those 
for the mid-thorax from group to group. Instead, fluences for the two locations 
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Figure 51. 
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Figure 57. 
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Figure 60. 
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oscillate around each other within the limits of their statistical uncertainty, 
which is of the order of 10 to 18 percent. At low energies there is little 
consistent difference between mid-head and mid-thorax for neutrons except at 


thermal, which translates into a similar difference in neutron-gammas (Figure 56) 
with the mid-thorax having the larger value. On the other hand the low energy 
gamma rays at the mid-thorax location are consistently higher than those for the 
mid-head for all energy groups less than an MeV. Since the simple phantom has no 
skull, skeletal shielding cannot be blamed for this phenomenon. It may be that it 
is due to the shielding of low energy photons by the wood box, which again puts a 
premium on the correct incident angular fluence description to calculate the 
correct transmission. 


The general trends for the simple phantom in room two (Figures 59 through 62) 
are similar to those described for room one. An exception to this is that, in the 
absence of the lead shield, the free-field gamma ray fluence is generally greater 
than that in the phantom at energies above one MeV. However, below that energy 
the build up from scattering of externally-incident gamma rays causes the 
transmitted fluence to exceed the free-field. 


In summary the detailed results of the calculations indicate that the 
dosimetric behavior of the simple and complex phantoms is qualitatively similar, 
though transmission magnitudes may differ between the two. Further, the results 
show that transmission into a subject even as smal] aS a rhesus monkey is affected 
by variations in incident spectra, particularly in the case of neutrons. Finally, 
the results show that both the neutron and the gamma ray spectra at locations 
within a phantom, are remarkably consistent, regardless of the incident spectra, 
and are different than those of the free-field, with the biggest differences 
arising in the low energy regions as a result of down-scatter buildup. 


3.2 Scalar Coupling 


Thus far this report has described the results for forward-adjoint fluence 
coupling via the VCS code, which is a complex process requiring a large computer. 
It would be convenient if it were possible to approximate this process using a 


Po ee ee ed 


scalar fluence coupling approach which is either not dependent on incident angle 
or marginally so. The mathematical expression for such a procedure is: 


Detector Response -{f @(E,2) K*(E)dQdE 
EQ 


in multigroup form this reduces to: 


Detector Response 5 sy DUE,» Wi)K*(E;) 
ToJj 
where K*(E;) is the KERMA-weighted adjoint fluence from an isotropic detector 
integrated over all angles at the leakage surface for leakage energy group i as 
described previously. In the DNA group structure, n = 58, which results in a 
convolution process manageable on a small computer or even a hand calculator. A 
better approximation might be to couple by hemispheres, as in: 


n 2d Ya 
Detector Response, _ =) L ® 2.) >: K*(E, 2,) 
i i EY Jk 


where @ and K* are the forward and adjoint fluence, respectively, integrated 
individually over the leakage surfaces which constitute those in the direction of 
the source (forward) or away from the source (backward). 


Both of the above scalar processes have been tried. The results are 
presented in Tables 11 through 15 in terms of transmission factors, along with 
those obtained from the VCS computation, which is considered to be the standard 
for comparison. The most obvious conclusion which can be drawn from a comparison 
of the results of the various coupling modes is that the scalar processes are not 
good approximations of the full VCS coupling scheme. This is true of both the 
full scalar (4) coupling and the hemispherical (27) coupling. The hemispherical 
coupling indicates the great importance of the fluence coming from the direction 
of the reactor. However, even that more complex approximation is unable to 
properly account for what is clearly a match between the direction of the forward 
fluence and the angular sensitivity of the monkey system, either simple or 
complex. This means that there is a significant decrease in both fluence and 
angular sensitivity to that fluence for directions other than a direct line 
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between the center of the reactor and the detector. If this is so, then the 
effect should be greater for the mid-thorax location than the mid-head or marrow 
since the mid-thorax angular sensitivity is likely to vary more with polar angle 
than those for the mid-head or for the marrow, which has components near the 
surface of the monkey. Tables 11 through 15 indicate that this is indeed the 
case. Further, it is indicated that neutrons are more affected by this phenomenon 
than gamma rays. This is to be expected since the gamma rays are more 
penetrating, thereby giving a more uniform sensitivity as a function of angle. 


Finally, it can be seen from the tables that even the neutron-induced gamma 
ray component exhibits an angular sensitivity, which is between those of the other 
two components in severity. Normally, the neutron-induced gamma ray component is 
not so angle sensitive, because it depends on the incident total neutron fluence 
rather than fluence energy. However, in this case it is possible that the wood 
box is stopping the lower energy neutrons from reaching the phantom; thus 
preventing the production of gamma rays in a geometrically advantageous position. 
As a result of this, the window in the box may play some part in defining the 
angular sensitivity of neutron-induced gamma ray response. 


In summary, it would seem that simple scalar coupling methods should be 
avoided in estimating transmitted dose in the AFRRI experimental configuration. 
It is also probable that attempting to find some other scalar coupling scheme 
sufficiently flexible to cover all situations would require more time than 
implementing a version of VCS optimized for smal} computer use. 
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Table 11. Results of forward-adjoint scalar fluence coupling ne Re 
(rads per KW-min) mid-head complex phantom HUN 


Scalar Coupling Mode Scalar/VCS 


Location Component 41 2nBack 2nForward 2,Total VCS (FSD) 4 On 


Mid-head, Complex Phantom, Room 1, 100 cm from. core-centerline, 2" Pb 


intervening. 
n 3.70 0.03 3.76 3.78 4.96 (.03) 0.75 0.76 
ney 0.27 0.01 0.26 0.27 0.32 (.07) 0.84 0.84 
Y 1.83 0.03 1.81 1.85 2.02 (.03) 0.91 0.92 


Mid-head, Complex Phantom, Room 2, 100 cm from core-centerline, bare room. 
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n 5.65 0.06 5.73 5.79 7.97 (.04) 0.71 0.73 
n-7 0.60 0.03 0.58 0.61 0.74 (.06) 0.81 0.82 ie 
Y 16.08 0.10 16.18 16.27 18.21 (.03) 0.88 0.89 » 
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Table 12. Results of forward-adjoint scalar fluence coupling 
(rads per kw-min) mid-thorax, complex phantom 


Scalar Coupling Mode Scalar/VCS 


Location Component 41 2iBack 2nForward 2nTotal VCS (FSD) 4, Qn 


Mid-thorax, Complex Phantom, Room 1, 100 cm from core-centerline, 2" Pb. 


n 2.68 0.02 2.92 2.94 4.62 (.03) 0.58 0.64 
n-y 0.50 0.02 0.51 0.53 0.69 (.05) 0.72 0.77 
Y 1.71 0.03 1,74 ali 2.11 (.03) 0.81 0.84 


Mid-thorax, Complex Phantom, Room 2, 100 cm from core-centerline, bare room. 


n 4.26 0.04 4.59 4.63 7.76 (.04) 0.55 0.60 

n-y 0.91 0.04 0.93 0.97 1.47 (.05) 0.62 0.66 

Y 15.06 0.09 15.52 15.61 19.24 (.03) 0.78 0.81 
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Table 13. Results of forward adjoint scalar fluence coupling 
(rads per kw-min) active marrow, complex phantom 
Scalar Coupling Mode Scalar /VCS 


Location Component 47 2nBack 2xForward 2nTotal VCS (FSD) 41 an 


Active marrow, Complex Phantom, Room 1, 100 cm from core-centerline, 2" Pb 


intervening. 
~ n 3.60 0.02 3.93 3.95 5.34 (.03) 0.67 0.74 
n-Y 0.31 0.01 0.32 0.33 0.37 (.05) 0.84 0.89 
Y 1.78 0.03 1.81 1.84 2.04 (.02) 0.87 0.89 


26S GI: 


Active marrow, Complex Phantom, Room 2, 100 cm from core-centerline, bare room. 


n 5.47 0.05 5.94 5.99 8.34 (.03) 0.66 0.72 
n-Y 0.65 0.03 0.67 0.70 0.94 (.05) 0.69 0.74 
Y 15.68 0.09 16.18 16.27 18.51 (.02) 0.85 0.88 
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Table 14. Results of forward-adjoint scalar fluence coupling 
(rads per kw-min) mid-head simple phantom 


Scalar Coupling Mode Scalar/VCS 


Location Component 4, 2nBack 2nForward 2ntotal VCS (FSD) 4n an 
Mid-head, Simple Phantom, Room, 1 100 cm from core-centerline, 2" Pb intervening. 


n 3.42 0.02 3.44 3.46 4.46 (.03) 0.77 0.78 
n- 0.35 0.01 0.35 0.36 0.45 (.06) 0.78 0.80 
Y 1.82 0.03 1.79 1.83 2.15 (.03) 0.85 0.85 


Mid-head, Simple Phantom, Room 2, 100 cm from core-centerline, bare room. 


n 5.27 0.06 5.25 5.31 7.13 (.04) 0.74 0.74 

n- 0.70 0.04 0.70 0.74 0.98 (.05) 0.71 0.76 

Y 16.00 0.10 15,96 16.05 19.17 (.03) 0.83 0.84 
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Table 15. Results of forward-adjoint scalar fluence coupling Shi det aie 
r (rads per kw-min) mid-thorax, simple phantom ott 
‘ SAE 
BS ee. 
: Sytat 
; Scalar Coupling Mode Scalar/VCS * 
i Location Component 41 2nBack 2nForward 2nTotal VCS (FSD) 4n 2a ~ 
cy, Mid-thorax, Simple Phantom, Room 1, 100 cm from core-centerline, 2" Pb = 
& Fe 
intervening. Sean 
a oy 
3 n 2.58 0.02 2.72 2.73 4.46 (.03) 0.58 0.61 ; 
WW 
4 n-y 0.50 0.02 0.50 0.52 0.69 (.05) 0.72 0.75 
+) Y 1.70 0.03 1.69 1.72 2.05 (.03) 0.83 0.84 
a Mid-thorax, Simple Phantom, Room 2, 100 cm from core-centerline, bare room. 
y n 4.11 0.04 4.30 4.34 7.51 (.04) 0.55 0.58 
i 
3 n-y 0.89 0.04 0.90 0.94 1.41 (.05) 0.63 0.67 
Y 14.97 0.09 15.03 15.12 18.95 (.03) 0.79 0.80 z... 
% e ; 
k: cee 
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: nee 
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SECTION 4 
VERIFICATION AND VALIDATION 


It is intended that, at some future time, KERMA and scalar fluence 
measurements be made by AFRRI in the simplified phantom at the locations 
calculated and presented in the previous section. In lieu of these it is still 
possible to make some comparisons with previously obtained experimental data as 
reported by arrri(}4) for the purposes of verifying and validating the 
calculational procedures described herein. 


AFRRI has made several measurements of the free-field gamma ray and total 
KERMA levels at several locations in the two exposure rooms with and without 
simple phantoms and other experimental apparatus. In particular, measurements 
were made using a small simple phantom of the same trunk height and diameter as 
the calculated simple phantom. The configuration of that phantom is shown in 
Figure 63, along with the outline of the calculational model. The two are very 
similar and are made of identical materials. 


The smal) phantom was irradiated in its configuration with chair and box in 
exposure room one, behind the two inches of lead, at a distance of 90 cm from core 
centerline. The small phantom was similarly exposed in the otherwise bare 
exposure room two at a distance of 125 cm from core centerline. The free field 
was measured at the mid-head and mid-thorax locations in the absence of the 
phantom but with the chair and box in place. The results of these measurements, 
along with the equivalent calculations are provided in Tables 16 and 17 for 
exposure rooms one and two, respectively. 


The calculations were performed by first computing the transmission to each 
location {in the chair and box in each room, but in the absence of the phantom. 
This transmission was normalized to the reported free field values, as measured. 
This normalization was then carried over to the calculation of KERMA within the 
simple phantom. The measured free-field values appear in the first column in each 
table. The two quantities actually measured are the gamma ray and total KERMA 
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values, each with a reported precision represented by a fractional standard 
deviation of 0.025. The neutron KERMA is a derived quantity, obtained by taking 


the difference of the other two quantities. Assuming that the reported KERMA are 
: the means of well-behaved normal distributions the variance of the neutron KERMA 
it ; ‘ ‘ 

e may be obtained according to the expression: (15) 

a4 2 a 2 2 

°neutron "total * gamma 

ql where o = MEAN*FSD. 

t 


For exposure room one the measured and derived neutron free-field KERMA all 
have approximately the same precision. However, in the case of room two the 
neutron KERMA exhibit a significant decline in precision. This is due to the fact 
that the gamma ray and total KERMA in room two are of the same order, whereas in 
room one the neutron and total KERMA were of the same order. 


The calculated phantom KERMA values appear in the third column of data in 
Table 16 and 17. The FSD values corresponding to each transmitted KERMA were 
obtained by again assuming a well behaved normal distribution of mean values as 
: represented by the measured free-field and the calculated transmission factor 
" given in column 6 and calculating their values using the expression: (15) 


2_. 2 2 2 2 2 2 
op % 95 + mM, oP) + m5 o 


where of is the oe as of ane product of the means, m, and m, which have ee 
respective variances, %, and G . As a result of these calculations, it can be eo 
seen that the calculated phantom KERMA for room one have a high level of e. e 
precision, as do the phantom gamma ray and total KERMA in room two. However, the 
precision of the calculated neutron KERMA in room two can be no better than that 
of the free-field to which it is normalized, and, thus, it is considerably worse 


than for the other components. 


‘at 


The calculated phantom KERMA may be compared with the measured values shown 
‘in column two of Tables 16 and 17. The agreement between al] the measured and 
calculated KERMA values, neutron and gamma ray, are very good for room one. This 
4 is gratifying because of the high precision attributed to all the measured 


109 


A rer ere FE eT RO aS ee og SION og ae) SS Ue eS CR et at ead he de 
SS ENS nr Ne eee 


see 


a Fy Sn Pon a a ie ee 
ASRS SSN RN GIR Ge SRG 


ges eetarte 
if 


= 
ves 


x 


ae he 
or ae 


a 


IRXe 


bate age Ate 


SAT pd | SAI 


quantities. For room two the agreement between calculated and measured gamma ray 
and total KERMA is similar to that achieved for room one. However, the measured 
and calculated mid-thorax neutron KERMA are considerably at odds. In spite of 
this apparent discrepancy a simple students t-test indicates that the two values 
are Statistically indistinguishable to 95% certainty, due to their large 
associated FSD values. Therefore, it is possible to say that the only really 
large difference between the calculated and measured phantom KERMA may be 
attributed to the inherent uncertainties in measured quantities, possibly 
compounded by the statiscal uncertainty of the calculations. 


In fact, as long as the measured gamma ray KERMA is a small fraction of the 
total, the uncertainty of the neutron KERMA will be of the same order as that of 
the total, which is likely to be very good. This is typified by the results for 
exposure room one, where the lead shield substantially reduces the gamma ray KERMA 
in favor of the neutron KERMA. On the other hand, the neutron KERMA in room two 
is small compared to the total and gamma ray KERMA values; in this case the 
variance of the neutron KERMA, being the sum of the variances of the two larger 
components, results in a standard deviation which is at least a large fraction of 
the neutron KERMA value. 


This discrepancy involving phantom neutron KERMA may be too large to be 
tolerated if calculated and measured quantities are to be used side by side for 
correlation with biological effects. Fortunately, other information is available 
which helps to shed some light on the origin of this discrepancy. If one looks at Bre ee 
the free-field spectra for the two rooms as shown previously in Figure 45, which ) 
have been validated by experimental means, one can see that the neutron spectrum 
above 100 KeV, i.e., that which contributes most to KERMA, is harder in room two 
than in room one. It should therefore be expected that the neutron transmission 
factor for identical locations within the phantom should increase going from room 
one to room two. The calculated values show such an increase, going from .54 in 
room one to .58 in room two for the mid-thorax location. However, for that same 
location the neutron transmission derived from the measurements shows a 
significant decrease going from room one to room two. Thus, a physical argument 
favors the calculated value. 


aaa 


be Perr 


Peres es, ee a ne el ee Le eee eee ee Sa ee 


The results of this effort to verify and validate the calculations described 
in this report are mixed due to the statistical uncertainties of both the 
experimentally measured and the calculated results but do seem to indicate a 
satisfactory reliability overall. This is particularly true under those 
circumstances which are most advantageous to the experimental methods. On the 
other hand, the results also suggest that, although extreme care was taken to 
obtain high precision in these experiments, the level of precision achieved was 
probably not sufficient to work reliably with room two in its bare configuration. 
Under such circumstances, the calculational approach may be found to be a 
particularly useful compliment to experimental methods. 
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APPENDIX 


This appendix contains tabulated values for adjoint leakage, free-field and 
transmitted fluence corresponding to those presented in graphical form in the 
body of the report. These fluence values are tabulated in the 37 neutron - 21 
gamma ray energy group format of the DNA Few Group Cross Section Library 
(DLC-31). The adjoint leakage fluence is tabulated in units analogous to KERMA 
factors, i.e., rads per unit fluence. Free-field and transmitted fluence values 
are tabulated as group fluence, i.e., total fluence per group, and ma: be 
convoluted with any set of KERMA by simple matrix multiplication to obtain total 
transmitted KERMA values. 
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ATTN: W. Murray 


National Institutes of Health 
ATTN: Library, Acquisition Unit 


National Science Foundation 
ATTN: Kin-Ping Wong 
ATTN: P. Haarlman 


Natl Heart, Lung & Blood Institute 
ATTN: W. Zukel 


National Library of Medicine 
ATTN: Library 


Office of Technology Assessment 
ATTN: M. Gough 


Office on Smoking & Health 
ATTN: J. Pinney 


US Senate 
Subcommittee of Nuc] Reg Comm for Envir & Public Works 
ATTN: J. Curtiss 


US House of Representatives 
Committee on Armed Services 
ATTN: Subcommittee on Mil Per & Comp 


US House of Representatives 
Committee on Veteran's Affairs 
ATTN: C. Graves 
ATTN: C. Moore 
ATTN: F. Stover 
ATTN: M. Fleming 
ATTN: R. Wilson 


US Public Health Service 
ATTN: Library 


US Senate 
Committee on Armed Services 
ATTN: J. McGovern 


US Senate 
Committee on Veterans Affairs 
: J. Steinberb 


ATTN: J. Susman 
ATTN: K. Burdick 
ATTN: T. Principi 
ATTN: V. Raymond 
ATTN: W. Brew 
US Senate 
Committee on Governmental Affairs 
ATTN: S. Ulm 


US Nuclear Regulatory Commission 
ATIN: R. Whipp for F. Arsenault 
ATTN: R. Whipp fpr R. Minogue 
ATTN: R. Whipp for W. Mills 


US Public Health Service Hospital 
ATTN: €. Nishimura 
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OTHER GOVERNMENT AGENCIES (Continued) 


Veterans Adm Wadsworth Hospital Center 
ATTN: T. Makinodan 


Veterans Admin Medical Center 
ATTN: K. Lee 


Veterans Admin Medical Center 
ATTN: 0. McGregor 


Veterans Admin Medical Center 
ATTN: C. Tessmer 


Veterans Administration 


: B. Polcari 
ATTN: D. Bosch 
ATTN: J. Donsbach 
ATTN: J. Smith 
ATTN: L. Hobson 
2 cy ATIN: D. Starbuck 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Adminstration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 
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OTHER GOVERNMENT AGENCIES (Continued) 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 
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ATTN: Director 
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ATTN: Director 
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Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Dfrector 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 
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ATTN: Director 
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OTHER GOVERNMENT AGENCIES (Continued) 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Oirector 


Veterans Administration 
ATIN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


Veterans Administration 
ATTN: Director 


The White House 


ATIN: Ofc of Policy Dev (DP) 


US House of Representatives 


Committee on Interstate & Foreign Commerce 
ATTN: Subcommittee on Health & Envir 


FOREIGN AGENCIES 


Canadian Embassy 
ATTN: Library 


EDF-RETN 1 
French Engr Bureau 
ATTN: Library 


Medical Enclave, Ansari Nagar 
ATIN: A. Taskar 


Japan-Hawaii Cancer Study 
ATTN: G. Glober 


French Engineering Bureau 
ATTN: M. Delpla 


McGill University 
ATTN: R. O'Seasohn 
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FOREIGN AGENCIES (Continued) 


Comitato Nazionale Per L'Energia Nucleare 
ATTN: Library 


Puerto Rico School of Medicine 
ATTN: Library 


United Kingdom Scientific Miss 
ATIN: Military Liaison for D. FakTey 
2 cy ATTN: Publications (for MRC) SQ 128 


OTHER 


Brookhaven National Laboratory 
ATTN: A. Brill, Medical Dept 
ATIN: E. Cronkite, Medical Dept 
ATTN: M. Bender, Medical Dept 
ATTN: Technical Library 
ATTN: V. Bond 


California Institute of Technology 
ATTN: £. Lewis 
ATIN: R. Christy 


i of Chicago 
N: P. Meier 


University of Colorado 
ATTN: Library 


Columbia University 
ATTN: A. Bloom 
ATTN: Library 


Columbia University 
ATTN: Div of Biostatistics 


Cornell University 
ATTN: W. Federer 


University of Drew 
ATIN: Library 


Medical College of Georgia 
ATTN: L. Stoddard 


Harvard School of Public Health 
ATTN: J. Bailor 
ATTN: Library 
ATIN: R. Reed 


Harvard School of Public Health 
ATTN: 8B. MacMahon 


Harvard University 
ATTN: W. Cochran 


University of Hawaii 
ATTN: Y. Matsumoto 


Indiana University 
ATTN: F. Putnam 


fowa State University 
ATTN: T. Brancroft 


Johns Hopkins University 
ATTN: A. Kimbal) 
ATTN: A. Lilienfield 
ATTN: R. Seltser 


Kansas University, Agri & Applied Sci 
ATTN: H. Fryer 


OTHER (Continued) 


Kingston Hospital 
ATTN: K. Johnson 


Memorial Hospital for Cancer and Allied Diseases 


ATTN: P. Lieberman 


Memorial Sloan-Kettering Cancer Center 
ATTN: J. Laughlin 
ATTN: P. Marks 


Merck, Sharp & Dohme Intl 
ATTN: A. Bearn 


University of Miami 
ATTN: P. Hodes 


University of Michigan Medical School 
ATTN: J. Neel 


University of Michigan 
ATTN: R. Cornel? 


University of Michigan 
ATTN: F. Moore 


Minnesota Dept of Health 
ATTN: 0. Lilienfield 


University of Minnesota 
ATTN: J. Bearman 
ATTN: L. Schuman 
ATTN: Library 


Natl Council on Radiation Protection & Measurements 


ATIN: W. Sinclair 


University of New Mexico 
ATTN: C. Key 
ATTN: R. Anderson 


New York University Medical Center 
ATTN: N. Nelson 


New York University 
ATTN: A. Upton 
ATTN: 8. Posternack 
ATTN: Library 
ATTN: M. Eisenbud 


University of North Carolina 
ATTN: 8B. Greenberg 
ATTN: Library for Dean 


Northwestern University 
ATTN: H. Cember 


Oak Ridge Assoc Universities 
ATTN: D. Lushbaugh 
ATTN: E. Tompkins 
ATIN: J. Totter 


University of Oaklahoma 
ATTN: P. Anderson 


University of Oregon 
ATTN: B. Pirofsky 


Pacific Northwest Laboratory 
ATTN: S. Marks 
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OTHER (Continued) 


Pennsylvania University Hospital 
ATTN: S. Baum 


University of Pennsylvania 
ATTN: P. Nowell] 


University of Pittsburgh 
ATTN: €. Radford 
ATTN: Library 


University of Pittsburgh 
ATTN: N. Wald 


Rochester University Medical Center 
ATTN: C. Odoroff 
ATTN: G. Casarett 


University of Rochester 
ATTN: L. Hempelmann 


Saint Francis Hospita? 
ATTN: R. Blaisdell 


Medical University of South Carolina 
ATTN: P. Liu 


University of Southern California 
ATTN: J. Birren 


Stanford University Medical Center 
ATTN: J. Brown 


Stanford University 
ATTN: L. Moses 


Stanford University Hospital 
ATTN: D. Dorfman 


Texas A & M University 
ATTN: R. Stone 


University of Texas at Austin 
ATTN: H. Sutton 


University of Texas 
ATTN: C. Cook 


University of Texas 
ATTN: R. Stallones 


University of Texas 
ATTN: W. Sutow 


University of Texas 
ATTN: G. Taylor 


University of Utah 
ATTN: Library 


University of Utah 
ATTN: C. Mays 
ATTN: E. Wrenn 
ATTN: L. Lyons 
ATTN: Library 


Vanderbilt University 
ATTN: R. Quinn 


University of Washington 
ATTN: D. Thompson 


OTHER (Continued) 


University of Washington 
ATTN: A. Motulsky 


University of Wisconsin 
ATTN: J. Crow 


Yale University School of Medicine 
ATTN: J. Meigs 
ATTN: Library 


DEPARTMENT OF ENERGY CONTRACTORS 


University of California 

Lawrence Livermore National Lab 
ATTN: L. Anspaugh 
ATTN: Tech Info Dept Library 
ATTN: Y. Ng 


Los Alamos National Laboratory 
ATTN: J. Dummer 
ATTN: Library 
ATTN: M/S634, T. Dowler 
ATTN: MS218, P. Whaten 


Oak Ridge National Laboratory 
ATIN: C. Richmond 
ATTN: G. Kerr 
ATTN: J. Auxier 


Oak Ridge National Laboratory 
ATTN: T. Jones 


Reynolds Electrical & Eng Co, Inc 
ATTN: J. Brady 
ATTN: LST 
2 cy ATTN: CIC 


Sandia National Laboratories 
ATTN: Div 1314, S. Durpee 


DEPARTMENT OF DEFENSE CONTRACTORS 


Advanced Res & Applications Corp 
ATTN: R. Armistead 


BDM Corp 
ATTN: J. Braddock 


Colorado State University 
ATTN: M. Zelle 


Energy Systems, Inc 
ATTN: T. Gates 


JAYCOR 
ATTN: A. Nelson 


Kaman Tempo 
ATTN: DASIAC 


LA Univ School of Medicine, Shreveport 
ATTN: Library 


National Academy of Sciences 
ATTN: National Materials Advisory Board 
ATTN: S. Jablon 
ATTN: S. McKee 
7 cy ATTN: C. Robinette 


DEPARTMENT OF DEFENSE CONTRACTORS (Continued) 


University of Nebraska 
ATTN: Library 


Ohio State University 
ATTN: Library 


Pacific-Sierra Research Corp 
ATIN: H. Brode, Chairman SAGE 


R&D Associates 
ATTN: C. Lee 
ATTN: P. Haas 


R&D Associates 
ATTN: A. Deverill 


Radiation Research Associates, Inc 
ATTN: N. Schaeffer 


Rand Corp 
ATTN: Library 
ATTN: P. Davis 


DEPARTMENT OF DEFENSE CONTRACTORS (Continued) 


Rand Corp 
ATTN: B. Bennett 


Science Applications Intl Corp 
ATIN: C. Thomas 
ATTN: J. Klemm 
ATTN: W. McRaney 
2 cy ATIN: J. Goetz 
5 cy ATTN: J. McGahan 


Science Applications, Inc 
ATTN: J. Striegal 


Science Applications, Inc 
2 cy ATTN: OD. Kaul 
2 cy ATIN: S. Egbert 
2 cy ATTN: J. Roberts 


Scientific Info Svcs, Inc 
ATTN: Library 


Varian Associates, Inc 
ATTN: E. Tochilin, Radiation Div C-063 
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